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In this paper, we present a few measurements done at the DO experiment at the Fermilab Tevatron 
Colhder. These measurements include the triple differential cross sections (d^a/dp^dy'^dy-''^*) of the 
photon and associated jet production, the photon and heavy flavour (6 and c) jet, and finally, study 
of the event with double parton scattering using 7+3 jet events. Each section below presents a brief 
description of those measurements and results. 



I. MEASUREMENT OF THE 
DIFFERENTIAL CROSS SECTION FOR THE 
PRODUCTION OF AN ISOLATED PHOTON 
WITH ASSOCIATED JET IN pp COLLISIONS 
AT v^=1.96 TEV 

The production of a photon with associated jets 
in the final state is a powerful probe of the dynam- 
ics of hard QCD interactions Different angular 
configurations between the photon and the jets can 
be used to extend inclusive photon production mea- 
surements and simultaneously test the underlying dy- 
namics of QCD hard-scattering subprocesses in differ- 
ent regions of parton momentum fraction x and large 
hard-scattering scales Q^. 

Here we present an analysis of photon plus jets pro- 
duction in pp collisions at a center-of-mass energy 
^/s =1.96 TeV in which the most-energetic (leading 
in pt) photon is produced centrally with a rapidity 
\y''\ < 1.0 [1]. The cross section as a function of 
photon transverse momentum p^ is measured differ- 
entially for four separate angular configurations of the 
highest Pt (leading) jet and the leading photon rapidi- 
ties. The leading jet is required to be in either the cen- 
tral (|yj<=*|<0.8) or forward (1.5 < lyj'^*] < 2.5) rapidity 
, intervals, with p!^^ > 15 GeV. The four angular config- 
urations studied are: central jets with y'^-y-'^* > and 
with ■ yj''* < 0, and forward jets with ■ y^'^^ > 
and with y'^ ■ y^°^ < 0. The total x and region 
covered by the measurement is 0.007 J; .t < 0.8 and 
900 < = (jPr^f < 1.6 X 10^ GeV^, extending 
the kinematic reach of previous photon plus jet mea- 
surements (see refs. inside 0). Ratios between the 
differential cross sections in the four studied angular 
configurations are also presented. The measurements 
are compared to the corresponding theoretical predic- 
tions. 

The data presented here correspond to an inte- 
grated luminosity of 1.01 ± 0.06 fb^^ collected using 
the DO detector [Si] at the Fermilab Tevatron pp col- 
lider operating at a center-of-mass energy ^/s =1.96 
TeV. 

Events containing at least one hadronic jet are se- 
lected. Jets are reconstructed using the DO Run II 
jet-finding algorithm with a cone of radius 0.7 and 
are required to satisfy quality criteria which suppress 
background from Icptons, photons, and detector noise 



effects. Jet energies are corrected to the particle level. 

In total, about 1.4 million candidate events are se- 
lected after application of all selection criteria. A cor- 
rection for the "7-l-jet" event purity V is then applied 
to account for the remaining background in the re- 
gion Onn > 0.7. The distribution of the ANN output 
for the simulated photon signal and dijet background 
samples are fitted to the data for each p^ bin using a 
maximum likelihood fit Q to obtain the fractions of 
signal and background components. The found puri- 
ties vary as 0.52 - 0.99 for 30 < < 400 GeV. 

The data are compared to next-to-leading order 
(NLO) QCD predictions obtained using JETPHOX 0] , 
with CTEQ6.5M PDF % and BFG fragmentation 
functions of partons to photons Q. The renormal- 
ization, factorization, and fragmentation scales {fJ-R, 
^p, and /i/) are set equal to p]^f{y*), where /(y*) = 
{[1 + exp(-2|y'^|)]/2}V2 ^^^^ y* = 0.h{y^ - yi-') [|. 

The prediction using the CTEQ6.5M PDF and 
BGF fragmentation sets does not describe the shape 
of the cross section over the whole measured range. In 
particular, the prediction is above the data for events 
with < 0.8 in the region > 100 GeV and be- 
low the data for jets produced in the 1.5 < < 2.5, 
y-f . yjot Q rapidity region for p^ < 50 GeV. Most of 
the data points in these and rapidity regions are 
(1-1.5) Satot outside of the CTEQ6.5M PDF set un- 
certainty range which is shown by the shaded region 
in the figure and calculated according to the prescrip- 
tion in [6] . Note that the data-to-theory ratios have a 
shape similar to those observed in the inclusive pho- 
ton cross sections measured by the UA2 [10], CDF [ll| 
and DO [l2| collaborations. 

The experimental systematic uncertainties are re- 
duced further by measuring the ratios between the 
differential cross sections D = d'^ a / dp^dy'^ dy''^^ in the 
different regions. Most of the systematic uncertainties 
related to the identification of central photons then 
cancel, and only systematic uncertainties related to 
the pp ^ ■^+}et + X event purities and the jet selec- 
tion efficiency (when measuring ratios between central 
and forward jet regions) remain. Measured ratios be- 
tween the differential cross sections in the different 
regions are presented in Fig. [51 The ratios also signifi- 
cantly reduces theoretical scale uncertainties [2] . The 
shapes of the measured ratios between the cross sec- 
tions in the different regions, in general, are qualita- 
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FIG. 1: The ratios of the measured triple-differential cross 
section, in each measured interval, to the NLO QCD pre- 
diction using JETPHOX 7] with the CTEQ6.5M PDF set 
and all three scales /ijj.F,/ = Prfiv*)- The solid vertical 
line on the points shows the statistical and pT-dependent 
systematic uncertainties added in quadrature, while the 
internal line shows the statistical uncertainty. The two 
dotted lines represent the effect of varying the theoretical 
scales by a factor of two. 

tively reproduced by the theory. A quantitative differ- 
ence, however, between theory and the measurement 
is observed for the ratios of the central jet regions to 
the forward 1.5 < lyj"*! <2.5, y^-y^°^ > region, even 
after the theoretical scale variation is taken into ac- 
count. The ratio between the two forward jet cross 
sections suggests a scale choice /ii?,F,/ — 2p^f{y*). 
However, the ratios of the central jet regions to the 
forward 1.5 < <2.5, y'^ ■y^'^^ < region suggest a 
theoretical scale closer to fiR,F,f — 0.5p'^f{y*). 

II. MEASUREMENT OF + b + X AND 
-y + c + X PRODUCTION CROSS SECTIONS IN 
pp COLLISIONS AT = 1-96 TEV 

Photons produced in association with heavy quarks 
Q (= c or b) in the final state of hadron-hadron inter- 
actions provide valuable information about the par- 
ton distributions of the initial state hadrons [T3lll3|. 
Such events are produced primarily through the QCD 
Compton-like scattering process gQ — s- jQ but also 
through quark-antiquark annihilation qq ^ jg 
jQQ. Consequently, j + Q + X production is sen- 
sitive to the b, c, and gluon (g) densities within the 
colliding hadrons, and can provide constraints on par- 
ton distribution functions that have substantial uncer- 
tainties [6] . The heavy quark and gluon content is an 



FIG. 2: The ratios between the differential cross sections in 
each y''^^ region. The solid vertical error bars correspond 
to the statistical and systematic uncertainties added in 
quadrature while the horizontal marks indicate the sta- 
tistical uncertainty. NLO QCD theoretical predictions 
for the ratios, estimated using jetphox, are shown for 
three different scales: iir,fj=Ho, 0.5^o, and 2/io, where 
MO ^Prfiv*)- 

important aspect of QCD dynamics and of the funda- 
mental structure of the proton. In particular, many 
searches for new physics, e-g. for certain Higgs bo- 
son production modes [3, [la, 113] , will benefit from 
the increased PDF precision with respect to the heavy 
quark and gluon content of the proton. 

This section presents the first measure- 
ments of the inclusive differential cross sections 
d3cr/(dpjd?/''dyj'=*) foi j + b + X and + c + X 
production in pp collisions, where y'^ and y^°^ are the 
photon and jet rapidities 1^. The results are based 
on an integrated luminosity of 1.02 ± 0.06 fb~^ 
collected with the DO detector at the Fermilab 
Tevatron Collider at = 1.96 TeV. The highest 
Pt (leading) photon and jet are required to have 
I?;''! < 1.0 and \y'°^\ < 0.8, and transverse momentum 
30 < < 150 GeV and p^^* > 15 GeV. This selection 
allows one to probe PDFs in the range of parton- 
momentum fractions (x) 0.01 ^ x ^ 0.3, and hard 
scatter scales of 9 xlO^ < Q'^ = {pl f < 2x10^ GeV 2. 
Differential cross sections are presented for two re- 
gions of kinematics, defined by y'^ ■ y'°^ > and 
y-f . yjct ^ These two regions provide greater sen- 
sitivity to the parton x because they probe different 
sets of xi and X2 intervals, as discussed in 0]. 

Photon selection criteria and methods for estimat- 
ing photon purities are similar to those used in Qj- 
At least one jet must be present in each event. Jets 
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FIG. 3: Distribution of observed events for Pup-jot after 
all selection criteria for the bin 50 < < 70 GeV. The 
distributions for the b, c, and light jet templates are shown 
normalized to their fitted fraction. Error bars on the tem- 
plates represent combined uncertainties from statistics of 
the MC and the fitted jet fiavor fractions, while the data 
contain just statistical uncertainties. Fits in the other 
bins are of similar quality. 



are reconstructed using the DO Run II algorithm [4| 
with a radius of 0.5. Light jets are suppressed using a 
dedicated artificial neural network (6- ANN) [l^ that 
exploits the longer lifetimes of heavy-flavor hadrons 
relative to their lighter counterparts. The leading jet 
is required to have a 6- ANN output value > 0.85. 

About 13,000 events remain in the data sample after 
applying all selection criteria. To estimate the photon 
purity, the 7- ANN distribution in data is fitted to a 
linear combination of templates for photons and jets 
obtained from simulated 7 -I- jet and dijet samples, 
respectively. An independent fit is performed in each 
bin, yielding photon purities between 51% and 93% 
for 30 < < 150 GeV . The fractional contributions 
of b and c jets are determined by fitting templates of 
-PHF-jct = -InDi^track to the data, where P^^^^^ is 
the probability that a track originates from the pri- 
mary vertex, based on the significance of the track's 
distance of closest approach to the primary vertex. 
Templates are used for the shape information of the 
PuF-jct distributions. For b and c jets these are ex- 
tracted from MC events while the light jet template is 
taken from a data sample enriched in light jets, which 
is corrected for contributions from b and c quarks. 

The result of a maximum likelihood fit, normalized 
to the number of events in data, is shown in Fig. [3] for 
50 < < 70 GeV , as an example. The estimated 
fractions of b and c jets in all bins vary between 
25-34% and 40-48%, respectively. The correspond- 
ing uncertainties range between 7-24%, dominated at 
higher by the limited data statistics. 

The measured differential cross sections are shown 
in Fig. 2] for 7 -|- b -I- X and 7 -I- c -I- X production as a 
function of for the jet and photon rapidity intervals 
in question. The cross sections fall by more than three 



FIG. 4: The 7 -|- b -j- X and 7 -|- c -|- X diflferential cross 
sections as a function of in the two regions y^'y''^^ > 
and y''y^°^ < 0. The uncertainties on the data points 
include statistical and systematic contributions added in 
quadrature. The NLO pQCD predictions using CTEQ6.6M 
PDFs are indicated by the dotted lines. 
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FIG. 5: The data-to-theory ratio of cross sections as a 
function of p^ for j + h + X and ^ + c + X in the re- 
gions y^'y'"^ > and y''y'°^ < 0. The uncertainties on 
the data include both statistical (inner line) and full un- 
certainties (entire error bar). Also shown are the uncer- 
tainties on the theoretical pQCD scales and the CTEQ6.6M 
PDFs. The scale uncertainties are shown as dotted lines 
and the PDF uncertainties by the shaded regions. The 
ratio of the standard CTEQ6.6M prediction to two models 
of intrinsic charm is also shown. 



orders of magnitude in the range 30 < pj, < 150 GeV . 
The statistical uncertainty on the results ranges from 
2% in the first bin to ~ 9% in the last bin, while 
the total systematic uncertainty varies between 15% 
and 28%. 

NLO pQCD predictions [2^ , with the renormaliza- 
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tion scale /i^;, factorization scale ^p, and fragmenta- 
tion scale /i/, all set to p7p, are compared to data in 
Fig. m The ratios of the measured to the predicted 
cross sections are shown in Fig. [S) 

The uncertainty from the choice of the scale is es- 
timated through a simultaneous variation of all three 
scales by a factor of two, i.e., to lJ.R,F,f = 0-5pl and 
2p'^. The predictions utilize CTEQ6.6M PDFs 0, and 
are corrected for effects of parton-to-hadron fragmen- 
tation. This correction for b (c) jets varies from 7.5% 
(3%) at 30 < pI < 40 GeV to 1% at 90 < < 
150 GeV . 

The pQCD prediction agrees with the measured 
cross sections for j + h + X production over the en- 
tire p^ range, and with ^ + c + X production for 
p^ < 70 GeV. For p'^ > 70 GeV, the measured 
7 -|- c -|- X cross section is higher than the prediction 
by about 1.6-2.2 standard deviations (including only 
the experimental uncertainties) with the difference in- 
creasing with growing p^. 

Parameterizations for two models containing 
intrinsic-charm (IC) have been included in CTEQ6.6 
[T^ , and their ratios to the standard CTEQ predictions 
are also shown in Fig. [SI Both non-perturbative mod- 
els predict an increased 'y+c + X cross section, and for 
the BHPS model [l3| it grows with p^. The observed 
difference may also be caused by an underestimated 
contribution from the g — > QQ splitting in the annihi- 
lation process that dominates for > 90 GeV [2l| . 



III. DOUBLE PARTON INTERACTIONS IN 7 
+ 3 JET EVENTS IN pp COLLISIONS AT 
= 1.96 TEV IN D0 

Many features of high energy inelastic hadron col- 
lisions depend directly on the parton structure of 
hadrons which is still not yet well understood at both 
the theoretical and experimental levels. Phenomeno- 
logically, the proton (or antiproton) may be viewed, in 
the first naive approximation, as an object composed 
of three light quarks (or anti-quarks). The study of 
this structure is founded mainly on the use of a sim- 
plified theoretical model which considers high energy 
inelastic scattering of nucleons as a process involving 
a single collision of one quark or gluon from one nu- 
cleon with one quark or gluon from the other nucleon. 
In this approach, additional "spectator" partons do 
not take part in the hard 2 — > 2 parton collision and 
form the so-called "underlying event" . 

Another, much less developed, approach is based 
on models in which there might be more than one 
hard interaction of parton pairs in one collision of nu- 
cleons. Since each incoming hadron is a composite 
object, consisting of many partons, such a probabil- 
ity should be non-zero. Models with multiple parton 
collisions have been considered in a few theoretical 
papers (see refs. inside [l^). It is obvious that the 



rate of events with double parton scattering (DPS) 
depends on how the partons are distributed within 
the nucleon. The form of the parton spatial distri- 
bution and possible correlations between partons is 
practically unknown. This information is hard to ob- 
tain within the present theoretical models based on 
perturbative QCD and makes a relevant measurement 
particularly important. 

Recently the information about a fraction of the 
double parton interactions has become very important 
due to growing Tevatron luminosity and upcoming 
LHC experiments. It opens the possibility to search 
for signals from new physics processes for which the 
DPS events may give a noticeable background. 

So far, there have been only four dedicated measure- 
ments studying double parton scattering: the AFS ex- 
periment in pp collisions at y/s = 63 GeV [2^ , UA2 in 
pp collisions at y/s — 630 GeV [24| and two by CDF 
in pp colhsions yG =1.8 TeV |25l.[26l|. 

In the present analysis, we use a sample of 
7-1-3 jets events collected by the DO experiment 
during Run Ila at the integrated luminosity of 1 fb~^ 
in pp collisions at ^/s = 1.96 TeV 221. We determine 
the fraction of the double parton interactions in a 
single pp collision and also the value of "effective 
cross section" cTcff. The latter allows, with the known 
cross sections tr^ and for two independent parton 
scatterings A and B, to calculate ctdps cross section: 

ctdps = (-L) 

CTcff 

Here the normalization factor is a parameter that 
can be directly related to the distance between par- 
tons in the nucleon. If the partons are uniformly 
distributed inside the nucleon (large CTcff) j (^dps "will 
be rather low and, conversely, it should increase for 
a highly concentrated parton spatial densities (small 
CTcff)- A better energy measurement of photons as 
compared with jets helps in separating the DP scat- 
terings and allows us to better fix the scale of the main 
hard interaction. 

In the current analysis, to extract CTcff we use 
a technique that operates only with quantities de- 
termined from data analysis and minimizes many 
theoretical assumptions [l^, that were used in 
some previous measurements. Towards this aim, we 
measure the number of DP 7 -I- 3 jets events and the 
number of 7 -|- 3 jets events with hard interactions 
occurring in two separate pp collisions. The latter 
class of events will be called double interactions (DI) . 
Assuming independent (uncorrelated) scatterings in 
the DP process [H [Hi, the DP and DI events should 
be kinematically identical and differ by just amount 
of "underlying" energy in one- and two- vertex events. 
The expression for the effective cross section CTcff can 
be written in the following form: 

A^Di A^c(l) ^ 

''^^ - iVD;2Ar(2) ^ ' 
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Here A'^oipp) is the number of DI(DP) events, Nc{l) 
and Nc{2) are the numbers of beam crossings with 1 
and 2 hard coUisions, and Uhard is the cross section of 
the hard pp interaction. We measure deff in three bins 
of the 2nd (ordered in pt) jet (p^^*^) which serves as 
a scale of the additional parton-parton interaction. 

We used data collected with the DO detector during 
Run Ila which corresponds an integrated luminosity 
of about 1.02 ± 0.06 fb~^. The data sets and 7+jets 
selection criteria, analogous to the previous pp ^ j + 
jet + X measurements [2, [l^, have been used in this 
analysis. 

The main background for the DP events is caused 
by the 7 + 8 jets events resulting from the single 
parton (SP) scatterings with hard gluon radiation in 
the initial or final state qg — > q^gg or qq — > gjgg. 
The fraction of the DP events is determined in this 
analysis using a set of variables sensitive to the 
kinematic configurations of the two independent scat- 
terings of parton pairs, specifically to the difference 
between between the px imbalances of two pairs in 
7 + 3 jets events: 




Two other variables, similar to ([3]), have also been 
used in the analysis [2^ for a cross-check. In the equa- 
tions above, pt(7,*) and pt(j, fc) are the px vectors 
of the total transverse momenta of the two two-body 
system, and Spxi^, i), Spxij, k) are the corresponding 
uncertainties. The pairs are constructed by grouping 
7 with 3 jets in three possible configurations. The 
configuration that gives the minimum S is selected 
for each of the 5- family variables. In most events S is 
minimized by pairing the photon with the leading jet, 
while the additional jets both come from dijet system 
or one of them is replaced by the radiation jet. 

The A5-family variables are allied to the S-family 
They are defined as an azimuthal angle between the 
Pt vectors of the pairs that give a minimum S value: 

A5 = A0(p?,'^'^*% p^^*-^^*^) (4) 

where p^'-''^*' and jf^*^'^'^*'' are the total pT balance 
vectors of the pairs (j + jeti) and (jetj +jetk) which 
give minimal S. 

The models for signal DP (and DI) events have 
been constructed directly from data mixing 7-I- >1 jet 
and dijet minimum bias events in the one- (and two- 
) vertex data !22i. Then the fractions of DP events 
have been found from analysis of shapes for AS dis- 
tributions in the signal events and data in the adja- 
cent p^^*^ bins for which we should expect different 
DP fractions [l^. Fractions of DI events have been 
found fitting the shapes for AS* distributions of sig- 
nal and background events to data. The obtained DP 



fractions are shown in Tables [H while the DI frac- 
tions vary as 0.189 ± 0.029 at 15 < p^^*^ < 20 GeV, 
0.137±0.027 at 20 < pi^*^ < 25 GeV and 0.094±0.025 
at 25 < j^^^^ < 30 GeV. 

TABLE L Fractions of DP events /dp found for the three 
pi^'^ intervals (GeV). 



jot2 

Ft 


15-20 


20 - 25 


25-30 


/dp 


0.466±0.041 


0.334±0.023 


0.235±0.027 



TABLE II: Effective cross section a^g (mb) found in the 
three pj^'^ intervals (GeV). 



Ceft 


15-20 


20 - 25 


25-30 


FT 


16.2 ±2.8 


13.8 ±3.1 


11.5 ±4.7 



The resulting values of a^g with total (systematic 
statistics) uncertainties are given in Table |TT] for 
the three p^^*^ bins. The total systematic uncertainty 
varies between 19% and 31% and is mainly caused by 
the determination of the DI and DP fractions as well 
as by the ratios of the DP/DI selection efficiencies. 

One can see that the obtained (Tcff values in 
different jet px bins agree with each other within 
their uncertainties. Using this fact and also that the 
uncertainties in different jet px bins have very small 
correlation, we can calculate the (Toff value averaged 
over the three jet px bins. It gives us 

CT^^"''- = 15.1 ± 1.9 mb. (5) 

It is worth mentioning that the obtained average 
value is in the range of those found in previous anal- 
ogous measurements d^l-jl^. This closeness may in- 
dicate a tendency for a stable behavior of tJcff with 
respect to the transverse momentum of the jet pro- 
duced in the second parton-parton interaction. 

IV. CONCLUSION 

The presented measurements of production cross 
sections for 7+jet, 7+heavy flavour jet and study of 
the events with double parton scatterings provide pre- 
cision tests of perturbative QCD as well as informa- 
tion on the fundamental structure of the nucleon; mo- 
mentum and spatial distributions of partons in the 
nucleon. In addition to these important aspects of 
QCD dynamics, this knowledge should be very help- 
ful for many searches for new physics since, for ex- 
ample, 7+jet and 7 + 6(c) jets are noticeable com- 
ponents for the background to many Higgs boson de- 
cay modes (e.g. h 77, bb) as well to extra dimen- 
sions with graviton (G ee, 77) and technicolor (e.g. 
u;rp —J. 77ro, — > jtt^ with — > 66, tt^ — + be, be) 
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[29|. Knowledge of parton matter structure and par- 
ton correlations being extremely valuable per se is 
also important for understanding multijet production 



characteristic for many physics final states, especially 
expected in many SUSY models [2]| . 
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